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ABSTRACT 

We investigate the prospects for discovering blazars at very high redshifts (z ;> 3 — 6) 
with the Fermi Gamma- Ray Space Telescope (Fermi), employing a model for the 
evolving gamma-ray luminosity function (GLF) of the blazar population. Our previous 
GLF model is used as a basis, which features luminosity-dependent density evolution 
implied from X-ray data on active galactic nuclei, as well as the blazar sequence 
paradigm for their spectral energy distribution, and which is consistent with EGRET 
and current Fermi observations of blazars. Here we augment the high-redshift evolution 
of this model by utilizing the luminosity function of quasars from the Sloan Digital 
Sky Survey (SDSS), which is well-constrained up to z ~ 5. We find that Fermi may 
discover a few blazars up to z ~ 6 in the entire sky during its 5-year survey. We further 
discuss how such high-redshift blazar candidates may be efficiently selected in future 
Fermi data. 
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1 INTRODUCTION 

High-energy gamma-ray astronomy has progressed drasti- 
cally after the launch of the Fermi Gamma-ray Space Tele- 
scope (Fermi). Its Large Area Telescope (LAT) is cur- 
rently observing the ent ire gamma-ray sky in the 0.02 - 
300 GeV energy range (|Atwood et all 120091 ) ■ The major- 
ity of extragalactic sources detected by Fermi are blazars, 
a subclass of active galactic nuclei (AGNs) dominated by 
broadband, nonthermal emission aris ing from relativist ic 
jets or iented close to our line of sight (|Abdo et al.ll2009al lbl. 
l2010al ld). The Fermi 11 month catalog reports the detec- 
tion of ~ 600 blazars up to z ~ 3 at Galactic latitudes 
\b\ > 10° (|Abdo et all 12010 J ). It is expected that Fermi 
will detect more than 100 blazars i n the near future (e.g. 
Narumoto fe Totanil |200d : iDermerj 120071 ; llnoue fc Totani 



iJHartman et al.l ll999T ) Q On the other hand, by analyzing 
the S wift/BAT b lazar sample in the redshift range z = 0.03 
- 4.0, lAiello et al.l (|2009T ) have recently suggested that the 
density evolution of luminous blazars peak at z — 4.3. Since 
Fermi has an order of magnitude better sensitivity and po- 
sitional accuracy at high G alactic latitudes compared to 
EGRET JAtwood et alJlboOSh . it is naturally expected that 
Fermi will see much deeper into the universe in gamma-rays. 
The purpose of this paper is to discuss expecta- 
tions for the highest redshift blazars that Fermi may 
discover. This requires reasonable knowledge of their 
gamma-ray luminosity function (GLF). The blazar 
GLF has been discu ssed from different perspectives in 
many papers so far (|Padovani et alj 19931; Stecker et al.l 

IChiang et al 



20091 ). 
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Before the Fermi era, the Energetic Gamma-Ray Exper- 
iment Telescope (EGRET) on board the Compton Gamma 
Ray Observatory detected ~50 blazars in total up to z ~ 3 
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1 iRomani et al.l |200J) and iRomanil J2006I ) have reported one 
EGRET blazar candidate at 2 ~ 5.48, but it has not been con- 
firmed so far by Fermi. 
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after IT09) and llnoue et all J20ld , hereafter ITM10) 
have recently developed a new model of the blazar GLF 
that accounts for the blazar spectral sequence, as well 
as luminosity-dependent density evolution implied from 
the AGN X-ray luminosity function (XLF), and which is 
consistent with the EGRET and current Fermi data (see 
Sj2]). However, the ITM10 GLF is uncertain for z > 3, 
because the current observed number of X-ray AGNs and 
gamma-ray blazars above z ~ 3 is insufficient to strongly 
constrain the model. On the other hand, the optical lumi- 
nosity function (OLF) of AGNs has been well constrained 
up to 2 ~ 5 iJRichards et al.l 1200a ). Therefore, here we 
consider a modified blazar GLF by combining the AGN 
XLF with the evolutionary constraints from the AGN OLF 
data. 

This paper is organized as follows. We introduce our 
updated GLF model in SJU In Sj3l we show predictions and 
candidate selection methods for high-redshift blazar obser- 
vations with Fermi. A discussion and summary is given in 
<J4] Throughout this paper, we adopt the standard cosmo- 
logical parameter set (h, Qm, £2a) = (0.7,0.3,0.7). 



72 unidentified gamma-ray sources at |6| > 10°. However, a 
very recent analysis of the Fermi extragalactic source distri- 
bution revealed that the surface density of blazars is ~0.12 
deg" 2 downtoF(> 100 MeV) = lxlO" 9 photons cm" 2 s"\ 
which would impl y ~ 5, 000 blazars in the entire sky 
iJAbdo et alJl2010bl ). Although this is ~ 8 times larger than 
the actual number of Fermi blazars, the difference may arise 
from corrections for the detection efficiency of the LAT. The 
discrepancy from our prediction is a factor of ~4 down to 
this sensitivity. One possible source of the discrepancy is 
differences in the adopted range of gamma-ray luminosi- 
ties, since our GLF is based on the EGRET sample that 
contained only ~ 10 low-luminosity BL Lac objects (see 
lAbdo et alj|2010bf ). Another reason for the discrepancy may 
lie in the correction factors for the detection efficiency, which 
could be a s large as 10 -1000 n ear the sensitivity limit, ac- 
cording to lAbdo et al.l (|2010bl ) . More detailed analyses of 
the Fermi data including information on the distributions of 
luminosity and redshift is necessary for clarification. If the 
true number of sources is greater than that in our model, it 
would only increase the prospects for finding high-redshift 
blazars with Fermi. 



2 HIGH REDSHIFT EVOLUTION OF THE 
BLAZAR GAMMA-RAY LUMINOSITY 
FUNCTION 

2.1 Blazar Gamma-ray Luminosity Function 

IT09 recently developed a model for the blazar GLF fea- 
turing the so-called luminosity-dependent density evolution 
(LDDE), based on the latest determination of the AGN XLF 
that clearly show the tendency of the more luminous objects 
to undergo their peak ac tivity periods at higher redshifts 
l|Ueda et alJ 12003 : iHa singcr et al. 2005). Another novel as- 
pect of IT09 was the account of the blazar sequence, which 
refers to the observed trend whereby the two characteristic 
frequencies at where the blazar spectral energy distribution 
(SED) peaks syste matically decr e ase as the b o lometric lu- 
minos i ty increases (IFossati et al. J 1997 L 1998 ^j_ Kubo et al.l 



\ 19981) ifPonato et all (l20qil ): lGhisellini et all fcoogl ). but i 
also lPadovani et al.l (2007)). The key parameters in the GLF 
model have been carefully determined to match the ob- 
served flux and redshift distribution of EGRET blazars by 
a likelihood analysis. Although the blazar sequence SED in 
IT09 was observationally constrained only up to the EGRET 
band of 30 GeV, this was extended to include the TeV band 
in ITM10, using published TeV blazar data. In this paper, we 
use the ITM10 GLF as our baseline model. Since the main 
modification of the blazar sequence SED in ITM10 was for 
J> 100 GeV, predictions for Fermi by ITM10 are similar to 
tnose by IT09 (see ITM10 for details). 

Including som e contributions from radio-quiet AGNs 
ijlnoue et al.l l2008l ). the extragalactic gamma-ray back- 
ground (EGRB) spectrum predicted by IT09 was found 
to be consi s tent w ith that actually observed by Fermi 
(JAbdo et al.l l2010dh . Furthermore, the expected number 
count for Fermi blazars was ~ 750 all-sky for the typ- 
ical 1-year survey sensitivit y limit of F(> 10 MeV) = 
3 x 10" 9 photons cm" 2 s" 1 JAtwood et alj|200g| ). This also 
seems to be consistent with the numbers reported in the 
Fermi 11-month AGN catalog, which are 596 blazars and 



2.2 Constraints on the High Redshift Evolution of 
Blazars 

The IT M10 GLF model is based on data fro m EGRET 
blazars (Hartman ct al. 1999) and X-ray AGNs (|Ueda et al.l 
120031 : lHasinger et alj|2005j ). the highest redshifts for both 
samples being z ~ 3. To address the evolution at z ^ 3, 
additional observational constraints are necessary. Optical 
surveys such as the Sloan Digital Sky Survey (SPSS ) have 
succe ssfully identified quasars up to z — 6.43 (|Fan et al.l 
l2003h. and the AGN O LF is well determined up to z ~ 
5 JRichards et all [JOOq . hereafter R06). Utilizing the R06 
OLF, below we modify the high-redshift evolution of our 
previous best-fit blazar GLF based on the AGN XLF of 
lUeda et all l|2003l . hereafter U03). 

The AGN XLF and OLF arc merged following the pro- 
cedures described in §5.4 of Richards ct al. (2005) and §3, 
4 in U03, which we briefly summarize. First, since the R06 
OLF concerns optically selected type-I AGNs only while the 
U03 XLF includes all types of AGNs, the latter is converted 
as follows. From Fig. 9 in U03, the fraction of optically se- 
lected type-I AGNs as a function of absorption column den- 
sity Nh can be characterized as 

C(JVh) 

( 0.96 (20.0 < log Ah < 20.5), 

= 1 8.504 -0.368 log A H (20.5 sC log A H < 23.0), (1) 
[ 0.04 (23.0 < log A H < 24.0). 

From Eqs. 8-10 in U03, the probability distribution of Ah 
for AGNs having luminosity Lx at redshift z is 



ii(Nh;Lx,z) 



2 - ttt^x, z ) ( 20 -° ^ lo S ^h < 20.5), 



l+e 

—M L ^ Z ) 

T^*P(Lx,z) 



(20.5 < log A?h < 23.0), (2) 

(23.0 ^logV H < 24.0)., 



where Lx is the 2-10 keV luminosity in units of erg/s, 
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Figu re 1. MAz = 2] luminosity function at each redshift as indicated in each panel. Circles indicate SDSS quasar OLF data 
(jRichards et al.l 120061) . Dashed, solid, and dotted curves show the U03 original XLF (p2 = 1.5), modified U03 XLF (p2 = 3.5), and 
R06 OLF, respectively. 



e =1.7, and 
i/)(Lx,z) - min{i/>n 

for which 



, max[0.47 - 0.1 (log L x - 44.0), 0])}, 

(3) 



1pma,x — 



1 + e 



3 + e 
Then the fraction of optically selected type-I AGNs is 



(4) 



/< 



opt. type — I 



r2i.Q 

(Lx, z) = / dN H C(Nh) x V {L X , z; JV H ). (5) 

J20.5 

Next, we convert Lx to Mi[z = 2], the if-corrected i- 
band AB magnitude at z — 2. First, we evaluate the 2 keV 
luminosity L2k c v in units of erg/s/Hz by assuming a photon 
index of T — 1.9. This is then extrapolated to L2500, the 
luminosity at 2500 A by solving the equations 



logL 2 



-0.137 log L2500 + 2.638, 

^2kcV 



log L25OO + Oox log 



?'2r,iKi 



(6) 
(7) 



where q ox is the spectral e nergy index between 2500 A 
to 2 keV l|Steffen et alj |2006). and z^kov and z/2500 are re- 
spectively 2 keV and 2500 A in units of Hz. Finally, we 
convert from lumin osity L2500 to magnitude Mi[z — 2] as 
|Oke fc Gunnlll983l ). 



M t [z = 2] = -2.5 log 



47rd 2 



48.60-2.51og(l + 2), (8) 



where d = 10 pc = 3.08 x 10 19 cm. Thus the U03 XLF can 
be combined with the R06 OLF. 

The power-law index p2 in the U03 XLF characterizes 
the density evolution as a function of z at high redshift. 
Although P2 — 1.5 was used in ITM10, here we change this 



to P2 = 3.5 in order to be consistent with the high-redshift 
R06 OLF data. Since P2 does not affect the low-redshift 
evolution, the XLF is not significantly altered below z ~ 3. 
Fig. [T] shows the AGN OLF at each redshift in terms of 
Mi[z — 2] in comparison with the R06 OLF data. With 
this AGN XLF, we reconstruct the blazar GLF with the 
ITM10 blazar sequence SED model, following the procedures 
of IT09. 



Note that there are some discrepancies between the 
OLF data and our OLF model at the brightest luminosi- 
ties for low redshifts. This might be due to insufficient accu- 
racy in our method of converting luminosities and luminosity 
functions. The correla tion of Eqs|6]and[7]is derived from a 
sample of 293 AGNs l|Steffen et ail 12009 ). and the conver- 
sion betwee n the XLF t o OLF is based on a sample of only 
247 objects (|Ueda et al.l l2003L In order to reliably apply our 
model to wider ranges of luminosities and redshifts, the sta- 
tistical uncertainties must be decreased through more data 
from deeper, larger X-ray and optical AGN surveys. 



The key parameters of our new blazar GLF model are 
(g, 71, k) - (4.42,1.07,1.92 x 10" 6 ), where g is the ratio 
between the bolometric jet luminosity and nuclear X-ray 
luminosity, 71 is the faint-end slope index of the GLF, and 
k is a normalization factor for the GLF (see Section. 3 of 
IT09 for details). Based on this model, we now discuss the 
prospects for observing high-redshift blazars with Fermi. 
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F lim (>100MeV) = 1x10 9 ph/cm^/s 

1TM10 

This paper 

T(0.1-l GeV)>2.3 

r>22.2 
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Figure 2. Expected cumulative redshift distribution of blazars 
detectable by Fermi above 100 MeV in the entire sky for a flux 
sensitivity limit F(> lOOMeV) = 1 X 10 _9 photons/cm 2 /s. Solid 
and dashed curves correspond to whether or not AGN OLF con- 
straints are taken into account. The dotted curve is for T > 2.3, 
and the dot-dashed curve is for T > 2.3 together with r > 22.2 



3 PREDICTIONS FOR HIGH-REDSHIFT 
FERMI BLAZARS 

3.1 Expected Redshift Distribution of Fermi 
Blazars 

Fig. [2] shows the expected cumulative redshift distribution 
of Fermi blazars above 100 MeV in the entire sky for a 
~5-year survey flux sensitivity limit of F(> lOOMeV) = 
1 x 10 -9 photons/cm 2 /s. The effects of intergalactic ab- 
sorption due to 77 interactions with diffuse background 
radiation fields are not included here, as th ey are not ex- 
pected to be important below 1 GeV (e.g. iGilmore et al.l 
120091 : llnoue et all I2010J 1. Neglecting high-redshift optical 
constraints, Fermi may be able to detect blazars up to 
z ~ 10. Taking into account their high-redshift evolution 
implied from the AGN OLF, we expect that Fermi will find 
some blazars up to z ~ 6 with the ~5-year survey sensitivity. 



3.2 Methods for selecting high-redshift blazar 
candidates 

It is expected that Fermi will find mo r e than 1,000 
blazars in total ifNarumoto fc Totanil 120061 : iDermerl 120071 : 
llnoue fc Totanjl20Qg| ). Although ~ 70% of the Fermi blazars 
have confirmed re dshifts, the other ~ 30% still do not 
(|Abdo et al.ll2010al L Furthermore, 10% of the high-la titude 
Fermi sources remain unidentified (|Abdo et alj|2010d ). This 
implies that the highest redshift blazars must be searched 
for among a large number of sources and distinguished from 
numerous unrelated ones. The following methods for their 
observational selection should be effective for this purpose. 

First, we select the sources whose fluxes are close to the 
Fermi sensitivity limit, as the high-redshift blazars are nat- 
urally expected to be faint. Furthermore, it should be useful 
to look for thos e showing time variabi lity, a characteristic 
trait of blazars IjUrrv fc Padovanilll995T ). although this may 
not always be practical due to the faint fluxes. 

Second, we choose the sources whose 0.1-1 GeV photon 



indices Y > 2.3 to pick out blazars with gamma-ray lumi- 
nosities L 1 > 10 47 erg/s at 100 MeV, since high-redshift 
objects that are detectable should be luminous, and their 
spectral indices are expected to be soft according t o the 
blazar sequence (|Fossati et al.lll99cJ ; lKubo et alJll998h . Al- 
though this criterion depends entirely on the assumed blazar 
sequence paradigm which is still a matter of debate, the 
above trend is already seen in the 3-month bright Fermi 
blazar catalog (|Ghisellini et alj|2009f ). 

Third, we search for the radio counterpart using deep 
radio survey catalogs, which should give us tighter con- 
straints on the source location and identification. According 
to the blazar sequence, the radio fluxes of blazars at z ~ 6 
would be ~ 20 mjy, fainter than the limiting sensitivity of 
radio-loud galaxy catalogs used for the current Fermi cat- 
alog such as the Combined Radio All-sky Targeted Eight 
GHz Survey C CRATES ; iHealev et alj|2007h. the Can didate 
Gamma- Ray Blazar Sur vey (CGRaBS; JHealev et alj|2008h . 
and the Roma-BZCAT (|Massaro et alj|2009h . More useful 
for our purposes are the Faint Ima ges of the Radio S ky at 
Twenty centimeters survey (FIRST; Iwhite et al.lll997t) an d 
the NRAO VLA Sky Survey fNVSS; ICondon et al.lll99J ). 
which reach down to 1 mjy at 1.4 GHz in the northern 
sky. Within the Fermi localization uncertainty of ~ 10' at 
|6| > 10°, the average radio source count is nominally ex- 
pected to be ~ 0.20 blazars, so finding the radio counterpart 
should be a crucial discriminant. 

Fourth, we select objects that are not detected in op- 
tical surveys such as the SDSS, in view of the unavoidable 
attenuation due to intergalactic HI for high-redshift sources. 
Since FIRST covers the SDSS survey area, a quarter of the 
entire sky, sources detected by the SDSS can be searched 
for their radio counterparts. Here we set the optical selec- 
tion criterion to be r > 22.2 in AB magnitude, the r-band 
limiting magnitude of SDSS. We caution that the SDSS sen- 
sitivities in the i and z band would be insufficient for detect- 
ing blazars at 2 > 3 even in the absence of intergalactic 
absorption. Thus, the SDSS data by itself cannot be used 
for identifying high-redshift blazars, but will nevertheless be 
valuable for rejecting low-redshift contaminants. 

The dotted and dot-dashed curves in Fig. [2] show cases 
where the above selection criteria have been imposed for the 
gamma-ray photon index and the optical magnitude. The 
expected total number of Fermi blazars would decrease to 
~ 60 sources, of which ~38% is expected to lie at z > 5, so 
this procedure should be effective in narrowing down high- 
redshift blazar candidates, albeit with some remaining low- 
redshift contamination. Finally, deeper optical and infrared 
follow-up observations of the candidate objects are war- 
ranted to accurately identify their respective counterparts 
and spectroscopically measure their actual redshifts. Note 
that for blazars in the southern sky, surveys such as that 
being undertaken by the South Pole T elescope at 1.4 an d 
2.0 mm with milli-Jansky sensitivities (IVieira et alj|2010l ). 
as well as the VST ATLAS and VISTA surveys in the 
optical and near-infrared, should prove to be likewise valu- 
able for their selection. Even more powerful would be future 
all-sky near-infrared surveys by satellite missions such as 



2 VST ATLAS: |http://astro.dur.ac.uk/Cosmology/vst atlas/ 

3 VISTA: http://www.vista.ac.uk/ 
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the Joint Astrophys ics Nascent Universe Satellite (JANUS; 
iRoming et al.l 120101, Wide-field Imaging Surveyor at High 
Redshift (WISH fl or the Wide-Field Infrared Su rvey Tele- 
scope (WFIRST- Tstern et al.llioTol : lGehrelsll2010h . 
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4 DISCUSSION AND SUMMARY 

In this paper, we studied the prospects for detecting the 
highest-redshift blazars with Fermi through a model for the 
blazar GLF, and how we can select such blazars among 
numerous other expected Fermi sources. The high-redshift 
evolution of our GLF model was constrained from the ob- 
served AGN OLF, which is well determined up to z ~ 5. 
Thus we found that Fermi may discover some blazars up to 
z ~ 6 down to the ~5-year survey flux sensitivity limit of 
F(> lOOMeV) = 1 x 10" 9 photons/cm 2 /s. If the true num- 
ber of sources turns out to be larger than the predictions of 
our basic model (IT08,ITM09), the number of high-redshift 
blazars should be correspondingly larger. 

High-redshift blazar candidates may be selected 
through the criteria that the source (1) is faint and shows 
time variability, (2) has a soft spectrum with F > 2.3, corre- 
sponding to luminous objects in the blazar sequence, (3) has 
a radio counterpart in deep survey catalogs such as FIRST, 
and (4) are not detected in optical surveys such as SDSS, 
which should remove low-redshift contaminants. Then these 
sources must be followed up by optical and infrared tele- 
scopes. 

GeV photons from high-redshift sources are expected 
to suffer intergalactic absorption due to electron-positron 
pair production interactions with the cosmic UV back- 
ground radiation. The resulting attenuation features in the 
spectra of such sources could therefore be a key probe of 
the poorly-understood UV background and its evolution, 
providing valuable insight into th e epoch o f cosmic reion- 
ization and early star formation (Oh 2001; Gil more et al.l 
12003 ; llnoue et alJl201CO . Discovering high-redshift gamma- 
ray blazars will constitute a crucial step towards this goal. 

The typical isotropic bolometric luminosity Lboi ex- 
pected from the blazar GLF is 10 48 erg/s. If luminous 
blazars with Lboi > 10 49 erg/s are found at z J> 6, the 
implied mass of the central supermassive black hole would 
be Mbh > 1O 9 M0, assuming Eddington-limited jet power 
and typical jet Doppler factors ~ 10. Given the difficulty of 
explaining the existence of such massive black holes at z ;> 6 
with Eddington-limited accretion initiated from stellar-mass 
seeds (jTanaka fe HaimarJl2009 ; IVolonterill2010l ). future con- 
straints on high-redshift luminous blazars may also provide 
important constraints on the growth history of supermassive 
black holes. 
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